Gamma-ray signatures of annihilation to charged leptons in dark matter substructure 



Matthew D. Kistler and Jennifer M. Siegal-Gaskins 
Center for Cosmology and Astro-Particle Physics and Department of Physics, Ohio State University, Columbus, Ohio 43210 

(Dated: March 17, 2010) 

Due to their higher concentrations and smaU internal velocities, Milky Way subhalos can be at 
least as important as the smooth halo in accounting for the GeV positron excess via dark matter 
annihilation. After showing how this can be achieved in various scenarios, including in Sommerfeld 
models, we demonstrate that, in this case, the diffuse inverse-Compton emission resulting from 
electrons and positrons produced in substructure leads to a nearly-isotropic signal close to the level 
of the isotropic GeV gamma-ray background seen by Fermi. Moreover, we show that IfESS cosmic- 
ray electron measurements can be used to constrain multi-TeV internal bremsstrahlung gamma rays 
arising from annihilation to charged leptons. 

PACS numbers: 95.35.-l-d, 95.85.Ry, 98.70.Sa, 98.70.-f 



I. INTRODUCTION 

The possibility of using dark matter f]\ to provide en- 
ergetic positrons and electrons [2 6] to explain the GeV 
positron excess Q seen by PAMELA Hj, the "ATIC 
bump" I, [l0| . and the less- anomalous e -I- e"*" spec- 
tra measured by the Fermi Gamma -ray Space Telescope 
(hereafter, Fermi) Jll| and HESS 0, El (in lieu of a 
nearby pulsar [l4jlq |) has sparked considerable interest 
(e.g., p7H28j ). Generating the required flux through 
annihilations in the smooth component of the Milky 
Way's dark matter halo requires both a much larger anni- 
hilation cross section than might be expected for a ther- 
mal rehc [29I and a large branching ratio to charged lep- 
tons. A variety of constraints (e.g., [30-36]) already ap- 
ply to a smooth halo explanation, and will likely tighten 
with new data from Fermi [37| . 

We examine the observational consequences of anni- 
hilations occurring within the dark matter substructure 
of the Milky Way. Substructure differs from the smooth 
halo in its spatial distribution, which is less centrally- 
concentrated [H, [13, and in its characteristic velocity 
dispersions, which are colder. We focus on the case of 
enhanced annihilation to charged leptons in "Sommer- 
feld" models, in which the annihilation cross section in- 
creases with decreasing relative velocity, and consider 
scenarios with a velocity- independent cross section. An- 
nihilations to charged particles necessari ly p roduce in- 
ternal bremsstrahlung (IB) gamma rays pOl - |44| . High- 
energy electrons and positrons will also produce gamma 
rays through various energy-loss processes. Beyond a few 
tens of kiloparsecs from the Galactic Center, the domi- 
nant loss process for electrons is the inverse-Compton 
(IC) upscattering of CMB photons. 

We calculate the expected high-latitude gamma-ray 
emission from IB and IC resulting from dark matter an- 
nihilation in substructure throughout the Milky Way. In 
particular, rather than considering only the serendipi- 
tous presence of a single, nearby, massive dark matter 
clump .45, 46] , we account for the collective emission from 
the entire subhalo population. Noting that electron- and 
gamma-ray-induced showers are difficult to distinguish in 



air Cherenkov telescopes, we discuss how constraints on 
TeV gamma-ray fluxes can be obtained from the cosmic- 
ray electron measurements made by HESS to limit anni- 
hilations to lepton pairs via their IB emission. 



II. ANNIHILATION IN SUBSTRUCTURE 

To describe their structural properties, we assume each 
individual dark matter subhalo to be described by a 
NFW density profile [13, 
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where is a scale radius and Ps a characteristic density. 
The differential luminosity (photons or particles per en- 
ergy per time) L of a subhalo, for an annihilation cross 
section (cru)o that is independent of velocity, is 
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where M is the subhalo mass, c = r^ir/rs is the concen- 
tration, /(c) = ln(l + c) — c/(l + c), and the particle 
physics-dependence of the annihilation rate is isolated in 
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with rriDM the dark matter particle mass and dN/dE the 
particle spectrum produced per annihilation. 

Numerical simulations find a relation between concen- 
tration and mass for subhalos that varies as a function of 
distance from the Galactic Center (3^,13^. This is a natu- 
ral consequence of tidal stripping, which more effectively 
removes mass from the outer regions of the subhalos while 
leaving the core relatively unscathed. Thus, for a given 
subhalo mass, subhalos located nearer to the Galactic 
Center will be more luminous than those at large radii. 
We adopt the modified Bullock et al. [i^ relation for 
low-mass halos, with radial dependence, from Ref. [i^ 
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where rgeid = 402 kpc is the radius where equal mass 
subhalos and field halos have the same concentration (see 
also [EO])- Using this relation with Ec^. ([2]), we can ap- 
proximate the differential luminosity of a subhalo of mass 
M at a radius r from the Galactic Center by 



L{M,r) = KC{M) 
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where we have defined 
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to describe the dependence of the annihilation rate on the 
structural properties of the subhalo. We note that the de- 
pendences on M and r are not formally separable, but 
are weak enough that Eq. ([5|) gives a reasonable approx- 
imation. The mass of Canes Venatici I [HH, assuming a 
NEW density profile and concentration c = 19.5, is used 
to normalize and A/q. 

We assume a power-law mass function for the subha- 
los [38, 39], dN/dM oc M"" with a = 1.9 and extrap- 
olate this relation to a minimum subhalo mass Mmin — 
10-6 Mq. Noting that {dN/dC) = {dN/dM){dM/dC), 
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FIG. 1: Angular distribution of the emission from dark matter 
annihilation {J'{tp) from Eq. (IIOI) ). Top panel: Galactic sub- 
structure assuming radial distributions from Aquarius (solid) 
and Via Lactea II (dashed), compared to the smooth halo 
(dotted) assuming NFW, Einasto (a = 0.17, r_2 = 20 kpc), 
and cored isothermal (rcoro = 5 kpc) density profiles (as la- 
beled). Bottom panel: The substructure emission assuming 
the Einasto radial distribution with Aquarius parameters, 
broken down into contributions from radial shells surround- 
ing the Galactic Center (with distances as labeled) . Peaks are 
due to shell boundaries. 



integration over the subhalo population yields 
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which contains the dependence of the annihilation rate 
on the structural properties and mass function of the 
subhalos, and is independent of position in the Galaxy. 

We model the subhalo number density (i.e., number 
of subhalos per volume) at a radius r from the Galactic 
Center with an Einasto profile [s^l 
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with a^^ij^ = 0.68 and r_2 — 199 kpc, as found by the 
Aquarius Project [s^. We assume a NFW profile for 
the smooth halo with = 20 kpc, r^n- = 255 kpc, and 
Afvir = 1-9 X 10^^ M0, and normalize the subhalo dis- 
tribution such that a fraction /subs — 0.15 of Afvir is 
bound in substructure. For comparison, we also con- 
sider a number density distribution as in Via Lactea II, 
"vL-iiW « (l + r/Rs)-^ ii, with Rs « 20 kpc. 

The intensity (/) of gamma-ray emission resulting from 
annihilation in substructure at an angle ip from the 
Galactic Center is 



/(^) = ^rQplJ(^), 
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in which vq — 8.5 kpc and J^{ip) is given by the linc-of- 
sight integral 
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where r — r(s,ip). This is shown in Fig. [T] for the Aquar- 
ius and VL-II subhalo radial distributions, along with the 
analogous term for three possible smooth halo density 
profiles. Here, we have chosen C^^^^^ such that the anni- 
hilation rate per volume is matched to the smooth halo 
at vq. Both subhalo models produce a nearly-isotropic 
angular signal, much less strongly-peaked towards the 
Galactic Center than that of the smooth halo. In the fol- 
lowing, we adopt the Aquarius distribution from Eq. ([S|) . 



III. SOMMERFELD ENHANCEMENT IN 
SUBSTRUCTURE AND THE SMOOTH HALO 

A challenge in attributing the positron excess to dark 
matter annihilation is the need for a much larger anni- 
hilation cross section than expected for a thermal relic. 
One way to accomplish this is by introducing a scalar or 
vector boson that mediates an intermediate-range force 
between dark matter particles and can dramatically en- 
hance the cross section for annihilation at low relative 
velocities (see Sommerfeld [H^ for electromagnetic scat- 
tering). For dark matter, this can be through Standard 
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Model gauge bosons 
particle (</>) 

We examine the possibility that annihilations in sub- 
structure, rather than in the smooth halo, are the dom- 
inant contributor to the measured local lepton flux and 
the astrophysical consequences that result. In this Sec- 
tion, we first show how this can arise for Sommerfeld 
models, then proceed to properties valid in general. The 
Sommerfeld enhancement can be expressed by the factor 
S{v) = {av)/{av)o, where v is the relative velocity of the 
dark matter particles and {(tv)q is the cross section in 
the absence of the additional force. These models can 
be parametrized by the coupling of the dark matter to 
the mediator (as) and the ratio of the dark matter and 



or a new mediator models without resonant behavior, as in 156 



mediator masses, m-DM/m^. For small velocities 17 1 
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The enhancement saturates at the velocity Wmin at which 
the de Broglie wavelength of the dark matter parti- 
cles becomes comparable to the range of the force, i.e. 
^/(jnuMv) ~ The maximum enhancement is 
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Resonances corresponding to bound states can result in 
significantly larger enhancements in such models, when 
"1dm/w0 ~ 1r? jas-, for integer n [19]. We account for 
this by solving the nonrelativistic Schrodinger equation 
as in [53|- For simplicity, we consider as — 10"^, 

and note that this choice qualitatively demonstrates the 
features of these models. 

We also consider a more general parametrization of 1 jv 
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FIG. 2: Subhalo and smooth halo contributions to the local 
annihilation rate as a function of rriDM/m,^ in models with 
resonant behavior, for as = 10~^. Top panel: Ratio of the 
annihilation rate in substructure to the smooth halo {solid, 
left axis) and total annihilation rate (dotted, right axis). Bot- 
tom panel: Local enhancement factor in the smooth halo. 
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where Vq and 5*0 are the relative velocity of smooth halo 
particles and enhancement factor at tq. Without allow- 
ing for resonant behavior, large values of mDM/m^ are 
required to produce large enhancements. 

In the smooth halo, the local annihilation rate per vol- 
ume, F^^^^^j^, depends on the dark matter density and 
velocity dispersion at the solar circle. We normalize the 
1-D velocity dispersion to the local rotation curve [STj . 
;irc/\/2 ~ 180 km s^^, and take ctid ^ v. Then 
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As subhalos are dynamically colder than the smooth halo, 
their smaller internal velocity dispersions can lead to 
preferential enhancement of the annihilation rate. We 
adopt a simple relation between subhalo mass and ve- 
locity dispersion, cr„ '-^ -^^sub (^^^ US); ^^'^ assign 
an enhancement factor S[v{Msuh)] to each subhalo of a 
given mass by approximating the relative velocity of the 
particles throughout the subhalo by CTi, ~ w. (The r de- 
pendence of Csub in Eq. Q leads to a mild variation that 
we neglect here.) We normalize the (Jv{Msuh) relation 
using again the mass and velocity dispersion of Canes 
Venatici I [5l|. For Sommerfeld models, Eq. ^ becomes 
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The value of C^^^^^ is then determined by Eq. ([7]) using 
C{M) as defined in Eq. ([T5|). The local annihilation rate 
per unit volume from subhalos is 
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The large local number density of subhalos {n^^^^{rQ) ~ 
2 X 10* kpc"'^) allows us to neglect the discreteness of sub- 
halos as sources and use the annihilation rate averaged 
over the subhalo population via C^^^^ . Fig. [2] compares 
the local annihilation rates in the smooth halo and in sub- 
structure for a Sommerfeld model with as = 10~^. For a 
significant region of the mBM/m^ parameter space, par- 
ticularly near resonances, the subhalo rate vastly exceeds 
the smooth halo rate (top panel, solid line). In Fig.[3l the 
local smooth halo and substructure annihilation rates for 
a model without resonant behavior are compared. The 
two contributions are comparable at Umin ~ 10 km/s, 
with subhalos dominating at lower Wmin- 

At high subhalo masses (> 1O''M0) this approach 
breaks down due to the low number density of these 
subhalos, so that the presence of a nearby clump would 
have to be accounted for (although having such within 
1 kpc is unlikely [38|]). Including subhalos of all masses 
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FIG. 3: Same as top panel of Fig. (2] but as a function of Vmin 
for a Sommerfeld model without resonant behavior. The total 
rate {dotted Ime, right axis) is calculated for 5*0 = 50 (both 
the smooth halo and subhalo rates scale linearly with Sq). 



TABLE I: The ratio of local annihilation rates in substructure 
and the smooth halo Fsubs/Fsmooth and the fraction of the halo 
mass bound in substructure /sub for various choices of mass 
function slope a and minimum subhalo mass Mmin, in the 
absence of Sommerfeld enhancement. The maximum subhalo 
mass is 10"'^'' M0 in all cases. 



Minin (Me) 



/sub 



1.9 
1.9 
1.9 
2.0 
2.0 
2.0 



10-^ 
lO"** 
10-12 
10-^ 
10"^ 
10-12 



0.148 
0.150 
0.153 
0.442 
0.500 
0.672 



-L subs 



0.0366 
0.0452 
0.0815 
0.300 
0.541 
3.48 



in the calculation of the local flux results in an over- 
estimate of of not more than 25% relative to as- 
suming a maximum mass of 10^ Mq without accounting 
for a Sommerfeld enhancement, and smaller in models 
with lower Wmin, a minor effect that we neglect. This 
can be understood from Eq. ([7]): in absence of a Som- 
merfeld enhancement, the total luminosity roughly scales 
cx / dMM°-^'^dN/dM cx M-°;°^ so that each decade 
of mass contributes nearly equally. (This has a much 
weaker dependence on the choice of lower cutoff than the 

The choices we have made in defining our substruc- 
ture model are conservative: using a steeper mass func- 
tion (e.g., a = 2) , a smaller minimum subhalo mass 
M^ain (e-g-, |63 - l64| ) a larger mass fraction in substruc- 
ture /sub, a higher normalization of the c{M) relation 
(as in Aquarius) jOr considering substructures-within- 
substructure [H, [s^ would all increase the subhalo 
annihilation rate relative to the smooth halo. The depen- 
dence of the ratio of the local annihilation rates in sub- 
structure and the smooth halo Fsubs/rsmooth on a, /sub, 
and Mmin is explored in Table HI assuming a velocity- 
independent cross section (i.e., no Sommerfeld enhance- 
ment). For a = 1.9 (2.0) the mass function is normalized 
so that 15% (50%) (e.g., [1^) of the host halo mass is in 
subhalos of lO^^ to 10^° M©. 

The range of Fsubs/rsmooth in Table U indicates that 
even in the absence of an enhancement arising primar- 
ily from Sommerfeld effects substructure can contribute 
significantly to the local annihilation rate, and hence no 
more than a modest enhancement in the substructure 
rate is typically required for this component to be lo- 
cally dominant (see also [goI . [61|). Current CMB con- 
straints [g^ [g^ restrict the saturation cross section in 
Sommerfeld models to within a factor of a few of the 
value required in order to explain the PAMELA and 
Fermi data by annihilation in the smooth halo, conse- 
quently only models in which the subhalo contribution 
is no more than a factor of a few smaller than that of 
the smooth halo without Sommerfeld enhancement could 



produce sufficient cosmic-ray fluxes without exceeding 
the allowed saturation cross section. Clearly, this con- 
dition can be satisfied for plausible subhalo parameters, 
which is valid even for dark matter models with other 
means of obtaining an enhanced annihilation cross sec- 
tion. We proceed assuming that annihilation to charged 
leptons in substructure dominates over the smooth halo 
and accounts for 100% of the anomalous fluxes. This 
condition is sufficient for the following purposes, which 
can be rescaled linearly with the local annihilation rate 
in substructure as desired. 



IV. INTERNAL BREMSSTRAHLUNG 

For pure leptonic final states, the only gamma-ray 
emission directly resulting from annihilations in sub- 
structure is internal bremsstrahlung (IB), e.g., XX ~^ 
£~^£~"f. We consider a few representative annihilation 
channels: direct annihilation into 2 /i and 2 r, and anni- 
hilation to 4 T through a new particle (xx ~^ 4'4' and 
each (j) — )■ 2t) , with cross sections required to explain the 
combined PAMELA/Fermi data [63, ill. We calculate 
the IB spectra for the two lepton cases as in Ref. [4]| . 

dE E TT ^ ' \ mj^ ^ ) 

where a; = 1 — E/mj^yi and Oem — 1/137. For ttidm ^ 
nifj^^r, E, this has the behavior dNiB/dE oc E~^ [40| . 
We similarly calculate the four lepton detailed 
in Ref. . For the 2 /i channel, we omit the gamma-ray 
contribution from muon decay, — e z/g t'/i 7, which is 
negligible for the scenario considered here. 

The IB gamma-ray intensity from annihilation in sub- 
structure aX Tp — 180° (the minimum of the dark matter 
signal) is shown in Fig. 2] for the above cases, along with 
gamma rays resulting directly from pionic tau decays in 
the 2r scenario (using DarkSUSY [70|). Considering a 
smooth halo model would result in signals smaller by a 
factor of '-^ 2 — 3 for all profiles at large angles. Although 
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FIG. 4: Isotropic gamma-ray signals resulting from dark matter annihilations in substructure (assuming an Aquarius number 
density profile). Left side: Inverse-Compton gamma-ray emission of the final state electron/positron population from annihi- 
lations at distances > 20 kpc from the Galactic Center {solid line) and including all radii (dotted). This can be compared to 
COMPTEL [9l[, EGRET [92|, and Fermi [93| diffuse gamma-ray data. Right side: Internal bremsstrahlung associated with 
the birth of charged leptons is shown for annihilation to two muons (ttidm = 1.6 TeV; dot-dashed), two taus (rriDM ~ 4 TeV; 
dark dashed), and two pairs (muM = 8 TeV; double-dot dashed). For the two tau case, we show the effect of including tau 
decays (light dashed). Cosmic-ray e~ -\- e^ measurements from HESS (triangles) [l^. [T^ act as upper limits on an isotropic 
gamma-ray flux (see text). 



we do not otherwise consider thern, models based on de- 
cays in the smooth halo (e.g., |7ll - l74| ) would give signals 
comparable to those shown in Fig. H) 

Directly measuring such a diffuse gamma-ray flux at 
TeV energies is presently challenging, in part due to the 
effective area of Fermi saturating with energy |33| • While 
ground-based air Chcrcnkov telescopes do not have this 
problem, the electromagnetic showers that they observe 
are quite similar for TeV electrons and gamma rays, mak- 
ing them difficult to separate. Based on observations of 
fields far from the Galactic plane, HESS has recently 
reported measurements of the e~ -\- e"*" spectrum into 
the TeV regime [H, [l^l- In principle, a nearly- isotropic 
~ TeV gamma-ray ffux from dark matter annihilation 
could result in an app arent feature in this spectrum. As 
noted in Refs. [H, [l^, there should be little contribution 
from extragalactic TeV gamma rays. 

In Fig. m we show the HESS -I- e+ spectrum, which 
can be regarded as a conservative upper limit on isotropic 
TeV gamma rays. The maximum gamma-ray fraction of 
this measurement is likely ^ 10% (although systematic 
uncertainties could result in as much as ^ 50%) p^ . 
It is likely that a dedicated analysis that accounts for 



the fields of view observed by HESS and determines a 
limit on the photon fraction in a given energy interval 
can strengthen these constraints. Improved understand- 
ing of the underlying astrophysical electron spectrum 
would also allow for tighter constraints from high-latitude 
emission, while avoiding uncertainties associated with 
TeV Galactic Center emission, which is highly profile- 
dependent [TSj and would not apply here if substructure 
is depleted near the center of the Galaxy. 

V. INVERSE-COMPTON GAMMA RAYS 

Absent a means of containing them, high-energy elec- 
trons resulting from annihilations will escape subhalos 
without difficulty. Far from the Galactic disk, the most 
important loss channel is inverse-Compton scattering on 
the CMB (we neglect the cosmic IR background, which 
has energy density a few percent that of the CM B 17611 ) . 
since the magnetic fields there should be small [73; UM 
and result in negligible synchrotron losses. To calculate 
the gamma-ray fiux, we must first find the equilibrium 
e~ -|- spectrum. We start from the diffusion- loss equa- 
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tion for a spectrum of relativistic electrons, ne{E) |79l.l8C 



— [b{E)n,iE)] 



Q{E), (18) 



where the diffusion coefficient, V, is assumed to be 
isotropic, Q{E) is the source term, and b{E) = bo E^ is 
the radiative loss term, with 60 — 0.3 x 10^^^ GeV^^ s^^ 
for the CMB (in the Thomson limit). For dark mat- 
ter, equilibrium can be assumed. In an isotropic system, 
there is no dependence upon V, since particle losses are 
compensated for by gains. At ~ 1 TeV, the electron 
cooling time is ~ 10^ yr, so that even if electrons prop- 
agate rectilinearly, they would only travel a distance of 
order the virial radius of the Milky Way. This is likely an 
overestimate, since their propagation should be affected 
by the halo magnetic field, although its structure and 
strength is uncertain. Considering the length scales rel- 
evant for electrons injected by annihilation in substruc- 
ture, we make the simplifying assumption that this halo 
magnetic field results in the IC losses occurring near the 
injection point (more care is needed for smooth halo sig- 
nals due to the steeper gradient in particle injection with 



radius (Sl , 
equation 



I). This reduces the problem to a continuity 



,{E)]^Q{E)., 



(19) 



which can be readily solved for a given injection spec- 
trum. While IB signals may vary greatly between annihi- 
lation channels, essentially all models that remain viable 
post-Fermi lead to nearly identical equilibrium electron 
spectra (up to uncertainties in the astrophysical spec- 
trum and propagation models) [g^. 111] . With generality, 
we consider dark matter with mDM = 2.35 TeV anni- 
hilating into two /i^ pairs (as in [i^). The calculation 
proceeds similarly to [8J] as 



(20) 



where Tq is the local annihilation rate per vol- 
ume (matched to that required to agree with 
PAMEL A/Fermi) and the resultant IC gamma-ray spec- 
trum per annihilation, dNic/dE, is calculated using the 
methods of Ref. [s^. In the inner Galaxy, synchrotron 
and IC losses on optical/IR photon backgrounds would 
result in a broad range of secondary photons [85|. We 
thus consider the signal resulting from annihilations oc- 
curring beyond 20 kpc from the Galactic Center, where 
IC on the CMB can be safely assumed to be the dom- 
inant energy-loss mechanism (based on modeling of the 
Galactic optical/IR photon field [sej). Using J^(180°) 
for r > 20 kpc yields the solid line in Fig. 21 which can 
be seen from the bottom panel of Fig. [T] to be nearly 
isotropic. Naively including radii interior to 20 kpc would 
result in the dotted line. 

For the scenario considered here, the IC spectrum hap- 
pens to peak at a similar energy to the pionic spectrum 



from cosmic-ray interactions [87|. We note that this IC 
signal retains less angular information concerning sub- 
structure (88l - [90| than direct gamma rays (such as IB). 
In comparing to isotropic gamma-ray data (9l| - [93j . we 
have made no attempt to account for other astrophysical 
contributions (see, e.g., [ol,!!!] for blazars). 

The velocity-dependence of the annihilation cross sec- 
tion in Sommerfeld models makes calculating the cos- 
mic signal in this scenario more complicated than in the 
standard picture [96j . This requires moving beyond the 



97, m) due to 



assumption of a constant boost (as in 
a dependence of the velocity dispersion and hence the 
cross section on halo mass. Also, the effects of baryons 
on dark halos vary with mass, since low-mass halos likely 
were never able to retain gas to form stars, and even 
in halos containing dwarf galaxies dark matter governs 
dynamics in the inner regions [9^. Although a detailed 
treatment of these matters is beyond our scope, a sim- 
ple estimate based on the total annihilation rate within 
the Milky Way halo (and scaling the amount of mass 
within substructure with host halo mass) suggests that 
this could be a factor of a few larger than the IC flux from 
substructure in the Milky Way halo only (as in Fig. S]) 
with a similar spectral shape. We note that the scenar- 
ios assumed in Refs. [13, [sl] are fundamentally different 
than considered here, since gamma-ray measurements of 
the inner Galaxy can strongly constrain annihilations in 
the smooth halo for the density profiles considered and 
a velocity dependence in the enhancement was not taken 
into account in those studies. 



VI. CONCLUSIONS 

If one takes the position that astrophysical resolutions 
of the positron excess are untenable, then to have a viable 
dark matter scenario requires invoking a relatively-large 
annihilation cross section. One possibility for achieving 
this goal is a velocity-dependent cross section due to the 
presence of a new medium-range force resulting in a Som- 
merfeld enhancement. In this study, we have examined 
the observational consequences of annihilation to charged 
particles in the context of a halo populated with dark 
matter substructure. 

In determining the dark matter annihilation signatures 
of Galactic substructure, it must be kept in mind that 
the microphysics (annihilation) occurs within kinemati- 
cally distinct subhalos so that their macroscopic distribu- 
tion sets many aspects of the problem. We have demon- 
strated that for a range of models, Sommerfeld-based or 
otherwise, annihilations in substructure, rather than in 
the smooth halo, are the dominant source of the locally- 
measured lepton flux. These substructure-dominated 
models imply associated IB and IC gamma-ray emission 
at high latitudes at a level accessible to (and possibly 
already in tension with) current observations. We have 
also argued that HESS TeV electron measurements can 
be regarded as limits on the isotropic TeV gamma rays 
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arising from IB. Importantly, these new prospective sig- 
nals can be tested with upcoming gamma-ray observa- 
tions by Fermi and air Cherenkov telescopes. 
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